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Previous  reports  suggest  that  the  choice  of  base  catalyst  affects  both  the  kinetics  and  structural  evolution 
of  SiO;  gelation  from  tetraethoxysilane  (TEOS).  Furthermore,  these  effects  have  been  attributed  to  the 
interaction  of  the  alkaline  cations  with  charged  silicate  fragments.  In  this  study,  NMR  analysis  of  alkali-metal 
nuclei  and  of  ®Si  directly  shows  that  during  the  condensation  of  silicon  alkoxides  in  basic  alcohol  solutions, 
cations  form  ion-pair  complexes  with  negatively  charged  silicate  groups.  It  is  determined  that  the  ion  pair 
is  formed  most  favorably  with  potassium.  In  seeking  a  predictive  model  of  complex  stability,  we  found 
it  most  successful  to  use  a  bond  orbital  electronegativity  equalization  model  developed  by  Bergman  and 
Hinze. 


Introduction 

Much  of  the  promise  of  sol-gel  technology  rests  on  the 
development  of  reaction  engineering  models.  Therefore, 
we  are  interested  in  a  detailed  understanding  of  the  ki¬ 
netics  and  mechanism  of  sol-gel  reactions.  One  particu¬ 
larly  important  aspect  is  how  the  catalyst  affects  both  the 
overall  kinetics  and  the  structural  evolution  of  reacting 
silicon  alkoxides.  For  instance,  the  use  of  NaOH  as  a 
catalyst  instead  of  NH4OH  results  in  a  slower  condensation 
and  a  more  complete  conversion  of  tetramethoxysilane  2 
The  gel  microstructure  and  properties  can  also  depend 
strongly  on  the  type  of  base  used.  For  example,  though 
base-catalyzed  silicon  alkoxide  gels  are  generally  not 
spinnable,3  we  have  observed  that  when  KOH  is  the  cat¬ 
alyst,  spinnable  gels  are  formed,  whereas  this  is  not  the 
case  when  the  catalyst  is  NaOH  or  NH4OH. 

McCormick  et  al.4  have  previously  used  NMR  spec¬ 
troscopy  to  show  that  alkali-metal  cations  participate  in 
the  condensation  of  silicate  anions  in  aqueous  solution. 
Such  a  mechanism  may  also  be  operative  for  sol-gel  con¬ 
densation  in  solution  containing  both  alcohol  and  water. 
In  basic  conditions,  the  hydrolyzed  alkoxide  deprotonates 
to  form  anionic  groups  much  like  those  found  in  aqueous 
silicate  solutions.  We  have  previously  found  that  the 
concept  of  ion  pairing  helps  to  account  for  apparent  ac¬ 
tivation  energies  of  silicon  alkoxide  gelation  reactions 
under  conditions  where  the  gelation  results  from  cluster- 
cluster  condensation.5"7 

In  this  study,  we  present  cation  NMR  evidence  of  ion- 
pair  complexation  in  base-catalyzed  sol-gel  systems,  com¬ 
prised  of  TEOS,  water,  and  ethanol.  The  cation  NMR 
lines  allow  us  to  deduce  whether  the  cation  is  coordinated 
to  ethanol,  water,  or  silicates.  Furthermore,  we  test  three 
theoretical  models  to  predict  the  stability  of  the  ion  pair 
that  is  detected  by  NMR  spectroscopy.  First,  we  consider 
a  simple  electrostatic  model  based  on  formal  charges. 
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(2)  Kelt*.  L.  W.;  Effinger,  N.  J.;  Melpoder,  S.  M.  J.  Non-Crytt.  Solids 
1988,  83.  363. 

(3)  Brinker,  C.  J.  J.  Non-Crytt.  Solids  1988,  100,  31. 

(4)  McCormick,  A.  V.;  Bell,  A  T.;  Radke,  C.  J.  In  Better  Ceramics 
through  Chemistry  III-,  Brisker,  C.  J.,  Ed.;  Material  Reeearch  Society 
Sprinf  Meeting,  Reno,  April  1988. 

(5)  Schaefer,  0.  W.  MRS  Bull.  1988, 8,  22. 

(6)  Aeeink,  R.  A.  In  Better  Ceramics  through  Chemistry,  Blinker,  C. 
J.,  Ed.;  Wiley;  New  York.  1984;  p  301. 
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through  Chemistry  IV;  Brisker,  C  J.,  Ed.;  Material  Reeearch  Society 
oprinj  Meeting,  San  Francieco,  April  1990. 
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Second,  we  use  a  model  developed  by  Livage  and  Henry.* 
These  authors  have  pioneered  the  application  of  electro¬ 
negativity  equalization  to  predict  sol-gel  precursor  re¬ 
activities  through  the  use  of  a  partial  charge  model  based 
on  Allred-Rochow  electronegativities  and  Sanderson's 
principle  of  electronegativity  equalization.  However,  this 
model  might  not  successfully  predict  the  current  results. 
Finally,  we  consider  a  model  proposed  by  Bergmann  and 
Hinze,1  which  is  based  on  bond  orbital  electronegativity 
equalization  (BOEE).  We  show  that  this  model  predicts 
complex  stability  in  accord  with  experiment. 

Theory 

NMR  Spectra  of  Alkali-Metal  Nuclei.  NMR  spec¬ 
troscopy  is  a  powerful  quantitative  tool  to  study  the 
chemistry  of  the  group  la  elements.9-3*  In  spite  of  high 
nuclear  quadrupole  moments,  solvated  alkali-metal  cations 
ordinarily  give  sharp  single  NMR  signals.  This  is  attrib¬ 
uted  to  the  high  symmetry  and  mobility  of  solvated  cat¬ 
ions.9-3*  However,  alkali-metal  ions  are  known  for  their 
ability  to  form  ionic  complexes.  The  nuclei  of  cations 
complexed  with  ionophores  (e.g.,  valinomycin,  crown  eth¬ 
ers,  cryptands)  or  strongly  binding  anions  can  relax  quickly 
as  the  cation  loses  its  symmetrical  environment  and  mo¬ 
bility.19'12  This  broadens  the  NMR  signal,  sometimes  even 
beyond  detection.  For  example,  Shporer  et  al.13  were  un¬ 
able  to  detect  the  signal  of  MNa  complexed  with  valino¬ 
mycin.  In  this  case,  the  “disappearance"  of  the  NMR 
signal  indicates  that  the  cation  has  complexed.  Alkali- 
metal  NMR  spectroscopy  has  also  been  used  successfully 
to  observe  cation  complexation  with  a  number  of  organic 
acids14-15  and  silicate  oligomers  in  basic  aqueous  solution.4 

Partial  Charge  Distribution  Calculation.  Livage  et 
al.*-16  have  employed  partial  charge  distribution  calcula- 


(8)  Livage,  J.;  Henry,  M.  Ia  Ultrastructurt  Processing  of  Advanced 
Ceramics;  Ulrich,  D.  R.,  Ed.;  Wiley:  New  York,  1988;  p  183. 
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(11)  Haynes,  D.  H.;  Preeeman.  E  Kowaiaky,  A  Biochemistry  1971, 
10,852. 

(12)  Dye,  J.  L;  Andrewe,  C.  W.;  Ceraeo,  J.  M.  J.  Phys.  Chem.  1978 
79  3076. 

(13)  Shporer,  M.;  Poupko,  E;  Lux,  Z.  Inorg.  Chem.  It72,  11,  2441. 
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(15)  Rechnitx,  G.  A;  Zamochnich,  S.  B.  J  Am.  Chem.  Soc.  1984, 86, 
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cessful  in  some  cases,  this  model  does  not  account  for 
changes  in  charge  distributions  in  response  to  fhangesm 
molecular  structure  and  bonding-  Moreover,  Allred  Ro- 
chow  electronegativities  do  not  include  such  important 
parameters  as  the  electronic  valence  state  and  orbital 
polarizability. More  elaborate  methods  to  calculate 
charge  distributions  have  been  proposed;  they  allow  a 
variation  of  the  electronegativity  with  orbital  hybridiza¬ 
tion.20-21  In  some  models,  a  partial  equalization  of  elec¬ 
tronegativity  is  used.22"24  BeTgmann  and  Hinze1  have 
proposed  a  model  to  obtain  charge  distributions  based  on 
the  orbital  electronegativities  calculated  by  Hinze  and 
Jafffi.19-28  They  considered  a  linear  relationship  between 
the  electronegativity  of  an  orbital  i  and  its  charge  q,  fa 
=  0,  -1  or  -2  for  an  empty,  singly,  or  doubly  occupied 
orbital,  respectively): 

Xi  *  b  +  2cqt  (1) 

where  b  *  (3/  -  A)/ 2,  c  *  (7  -  A)/2,  and  /  and  A  are  the 
ionization  potential  and  electron  affinity  of  the  considered 
orbital-  For  each  two-electron  bond  A-B  formed  by  the 
orbital  i  of  the  atom  A  and  the  orbital  j  of  the  atom  B,  the 
equalization  of  electronegativity  is  given  by 

X,(<7;  +  A <?,)  *  Xjfa  +  H)  (2) 

where  Aq,  and  Aqt  are  the  charge  shifts  on  the  orbitals  i 
and  j. 

To  conserve  charge 


For  n  bonds,  a  system  of  n  linear  equations  is  obtained 
and  solved  numerically  for  the  Aq/g.  The  total  net  charge 
on  the  originally  neutral  atom  A  will  be 

Q( A)  -  £Aq,  (4) 
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The  model  can  be  improved  by  considering  a  linear 
dependence  between  the  parameter  b  and  the  “rest  charge" 

rr 

b  =  f>„  +  (5) 


where 


r,(A)  =  Q( A)  -  Aq, 


Table  I  gives  the  values  of  do,  6lf  and  c  used  in  our 
calculations. 

This  method  accounts  very  well  for  the  inductive  and 
polarizing  effects  of  chemical  groups.  It  also  distinguishes 
between  different  molecular  structures  or  isomers. 

Correlation  of  Partial  Charge*  with  Bond  Energy 
and  Bond  Length.  In  the  case  of  reactions  involving  polar 
interactions,  partial  charges  allow  prediction  of  the  re¬ 
activity  of  atoms  or  sites  in  relation  with  their  nucleophile 
or  electrophile  character.  Partial  charge  calculations  can 
also  be  used  to  estimate  the  Coulombic  contribution  to  a 
bond  energy,  which  is  given  by 

£c=P  (7, 

4»t0r 

where  Q*  and  Q~  are  the  partial  charges  on  respectively 
the  positively  and  negatively  polarized  atom  in  the  bond, 
<o  is  the  vacuum  permittivity,  and  r  is  the  bond  length. 
According  to  Derflinger  and  Polanski,28  the  bond  length 
is  approximately  determined  by  the  shortening  of  the 
purely  covalent  bond  length  by  the  charge  shift  Shomaker 
and  Stephenson29  have  shown  that  the  bond  length  d & 
can  be  estimated  with  the  relation 

<*AB  =  rA  +  rB  ~  *|XA  "  XbI  (8) 

where  rA  and  rB  are  the  covalent  radii  of  the  atoms,  xa  and 
Xb  are  their  electronegativities,  and  k  is  a  proportionality 
constant  between  the  bond  shortening  by  electron  transfer 
and  the  electronegativity  difference  of  the  atoms.  Berg- 
mann  and  Hinze1  determined  that  k  does  not  depend  on 
the  bond  type  and  it  is  a  constant  equal  to  1.9  pm/ V  (when 
the  electronegativities  are  in  volts).  The  above  formula 
has  successfully  predicted  bond  lengths  for  more  than  120 
compounds  within  a  standard  deviation  of  4  pm.  The 
covalent  radii  given  by  Sanderson80  have  been  used  in  the 
eq  8  for  our  calculations. 
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Figure  2.  “Si  NMR  spectrum  for  the  system  (in  molar  ratios) 
lTEOS/4EtOH/  lHjO/O.lNaOH  after  4  h  of  reaction:  (a)  TEOS 
signal;  (b)  (Et0)3-Si-O-Si-(0Et)3. 

Experimental  Section 

Tetraethoxysilane  (TEOS),  ethanol,  wateT,  and  a  base  catalyst, 
either  NaOH  or  RbOH,  were  mixed  with  various  molar  ratios  in 
polypropylene  vessels  at  20  °C.  After  the  sample  was  transferred 
to  a  10-mm  glass  NMR  tube,  “Na  NMR  spectra  were  acquired 
at  132.283  MHz  with  100  pulses  of  31.5  ms  and  a  delay  of  1  s 
between  pulses.  Chemical  shifts  were  referred  to  NaOH  50  wt 
%  in  water. 

^Rb  NMR  spectra  were  acquired  at  163.634  MHz,  with  1000 
pulses  of  28  ms  and  a  delay  of  0.1  s  between  pulses.  Chemical  shifts 
were  referred  to  RbOH  50  wt  %  in  water. 

“Si  NMR  spectra  were  recorded  at  99.348  MHz  with  100  scans 
of  25  ms  and  a  pulse  delay  of  3  s.  Peaks  were  referred  to  tetra- 
methylsilane  (TMS).  Cr(acac)3  was  used  as  a  paramagnetic  re¬ 
laxing  agent  to  reduce  the  spin-lattice  relaxation  time  of  “Si  and, 
thus,  to  ensure  quantitative  spectra. 

Results  and  Discussion 

NMR  Evidence  of  Ion  Pairing.  Figure  1  shows  “Na 
spectra  at  various  points  in  the  sample  preparation.  Before 
TEOS  addition,  the  “Na  line  is  narrow,  indicating  that 
the  cations  are  fully  solvated  and  are  moving  freely  in  the 
solution  (top  spectrum).  The  line  width  indicates  that  it 
is  solvated  by  water  rather  than  ethanol.38  Immediately 
upon  addition  of  TEOS,  the  sodium  signal  both  diminishes 
and  broadens  drastically,  either  because  of  the  introduction 
of  asymmetry  in  the  cation  environment  or  because  of  a 
lowered  mobility  resulting  from  complexation.  As  the 
reaction  proceeds,  the  cation  resonance  virtually  disappears 
(bottom  spectrum).  The  corresponding  “Si  NMR  spec¬ 
trum  (Figure  2)  shows  no  hydrolyzed  TEOS.  Thus,  as 
TEOS  is  hydrolyzed,  it  quickly  polymerizes  to  dense 
particles  such  as  those  detected  by  X-ray  scattering*  and 
electron  microscopy.31  “Si  signals  from  these  particles 
do  not  readily  appear  in  these  high-resolution  spectra 
because  (i)  the  fully  condensed  (04)  and  surface  (Q3)  Si 
sites  are  broad  due  to  chemical  shift  distributions,  (ii)  they 
would  relax  slowly  because  of  limited  Cr3+  access,33  and 
(iii)  they  are  partly  hidden  by  the  broad  signal  from  the 
glass  NMR  tube.  We  deduce  that  the  sodium  cations 
become  complexed  with  negatively  charged  groups  (de- 
protonated  Q3  silanols)  on  the  colloidal  particles.  The 
absence  of  any  narrow  “Na  signal  proves  that  no  Na+ 
cations  are  free  in  solution,  although  substantial  amounts 
of  TEOS  remain  (even  well  after  gelation). 

RbOH  catalyzed  systems  exhibit  the  same  kind  of  be¬ 
havior  (Figure  3).  Incidentally,  we  note  that  as  the  con¬ 
densation  reaction  nears  completion,  Si-O"— Rb4-  bonds  are 


(31)  Baiky,  R.  K.;  Nagaa*,  T.;  Brobarg,  S.  M.;  Macartney,  M.  L., 
Univanity  of  Minnesota,  personal  communication,  1990. 

(32)  Engelhanit,  G.;  Radeglia,  R.  Chtm.  Pf »y*.  Lett.  1994, 108,  271. 
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Figure  3.  87Rb  NMR  spectra  for  the  systems  (in  molar  ratios) 
5EtOH/ lH2O/0.1RbOH  (top  spectrum)  and  lTEOS/4EtOH/ 
lHjO/O.IRbOH  after  5  min  of  reaction  (bottom  spectrum). 
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Figure  4.  *7Rb  NMR  spectra  for  the  systems  (in  molar  ratioe): 
(a)  5EtOH/ lHjO/O.OIRbOH;  (b)  lTEOS/4EtOH/lH,0/ 
O.OIRbOH  after  5  min  of  reaction;  (c)  same  as  4b  but  with  LiCl 
at  saturation;  (d)  same  as  4b  but  with  NaCl  at  saturation;  (e)  same 
as  4b  but  with  NH4CI  at  saturation;  (f)  same  as  4b  but  with  KC1 
at  saturation;  (g)  same  as  4b  but  with  CsCl  at  saturation. 

replaced  by  Si-O-Si  bonds.  Evidently  the  Si-O-Si  bond 
formation  is  more  favored,  but  neither  water  nor  ethanol 
can  compete  to  form  Si-OH  or  Si-OEt  bonds. 

Summarizing,  base  cations  quickly  form  stable  Si- 
O'— M+  bonds,  even  at  low  cation  concentration.  These 
bonds  are  favored  relative  to  Si-OH  and  Si-OEt,  but  they 
are  unfavored  relative  to  Si-O-Si  bonds. 

Relative  Equilibrium  of  the  Cation  Binding.  The 
relative  equilibrium  of  various  ion  pairs  was  investigated 
by  adding  alkali-metal  chloride  salts  to  RbOH-catalyzed 
systems.  The  corresponding  ^b  NMR  spectra  are  shown 
in  Figure  4.  Parts  a  and  b  of  Figure  4  are  the  spectra  of 
the  system  with  and  without  TEOS;  as  previously  shown, 
the  Rb  peak  disappears  when  TEOS  is  added.  If  LiCl  is 
added,  no  change  occurs  (Figure  4c).  If  NaCl  is  added,  a 
small  amount  of  Rb4  is  released  into  solution  (Figure  4d). 
With  NH4CI  a  larger  amount  is  released  (Figure  4e).  If 
KC1  is  added,  moat  of  Rb4,  is  released  (Figure  40-  Finally 
if  CsCl  is  added,  only  a  small  amount  of  Rb4,  seems  to  be 
displaced  (Figure  4g).  We  deduce  from  Figure  4f  that  most 
of  the  Si-O*— Rb4,  bonds  have  been  replaced  by  Si-O'-K4' 
bonds. 
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Cation  Binding  with  Sol-Gel  Silicates 


These  spectra  suggest  the  following  order  of  efficacy  of 
the  cation  at  competing  for  charged  surface  silicate  sites: 

Li+  <  Na+  <  Cs+  =  NH4+  <  Rb+  <  K+  (9) 

Rationalizing  Ion  Pair  Stabilities.  K+  has  been  ob¬ 
served  to  produce  the  fastest  condensation  of  monomeric 
silicic  acid37  and  produces  the  fastest  base-catalyzed  TEOS 
gelations,  yielding  a  spinnable  gel.  We  propose  that  cation 
complexation  may  cause  these  results.  To  explain  the 
relative  equilibrium  of  cation-silicate  pairs,  we  next  pursue 
a  theoretical  investigation. 

The  ion-pair  reaction  can  be  written 

M+~OH2  +  -OSi  =  M+-'OSi+  H20  (10) 

where  M+  is  Li+,  Na+,  K+,  Rb+,  or  Cs*.  The  reaction 
represents  the  competition  between  a  surface  silicate  group 
and  a  water  of  hydration  for  the  position  next  to  the  cation. 
Note  that  we  are  ignoring  other  waters  of  hydration,  as¬ 
suming  that  they  will  not  affect  the  reaction  enthalpy 
change  (for  the  entropic  effect,  see  below).  The  physical 
implications  of  this  simplistic  model  are  considered  shortly. 
We  are  also  assuming  that  no  water  of  hydration  comes 
between  the  cation  and  the  silicate  fragment.  This  has 
been  suggested  by  McCormick  et  al.4  The  equilibrium 
constant  for  ion  pairing  (ip)  is 
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Figure  5.  Structure  of  the  silicate  anion  SijOjo*'  used  in  our 
charge  distribution  calculations.  The  oxygen  atoms  involved  in 
the  complexation  with  alkali-metal  cations  are  filled  in  black. 


Table  II.  Partial  Charge  Distribution  for  the  Ion-Pair 
Complexes  from  the  Cation  M  (M  *  Li,  Na,  K,  Rb,  Cs)  and 
the  Silicate  Fragment  of  Figure  S* 


M 

Q(M) 

Q(0)‘ 

Q(0)c 

(?(Si) 

Li 

747 

-550 

-408 

415 

Na 

768 

-558 

-411 

407 

K 

842 

-587 

-422 

379 

Rb 

859 

-594 

-425 

373 

Cs 

88 8 

-605 

-429 

362 

Kjp  =  e-^'IKT  (11) 

The  free  energy  of  reaction  is  given  by 

AG,  =  A Ht-  TAS,  =  AH,(M+~-OSi)  +  AH,(H20)  - 

AHf(M+-H20)  -  A/ffCOSi)  -  TAS,  (12) 

We  assume  that  the  entropy  change  of  this  reaction  is 
similar  for  different  alkali-metal  cations  since  the  reaction 
entropy  will  be  dominated  by  changes  in  the  arrangement 
in  water  molecules  of  hydration.  The  nature  of  this  change 
is  the  replacement  of  hydrating  water  by  the  silicate 
fragment,  and  we  assume  that  the  remaining  hydrating 
water  molecules  are  not  affected  by  the  ion-pair  com¬ 
plexation.  Therefore,  only  AHf(M+— 'OSi)  and 
-OH2)  change  in  response  to  a  change  of  the  cation.  These 
two  terms  are  directly  related  to  £(M+~OH2)  and  £- 
(M+— "OSi),  the  energies  of  the  bonds  M+—OH2  and 
M+— ~OSi.  Therefore 

A#f(M+~-OSi)  -  AHm*~H20)  *  E( M+~OH2)  - 

E(M+~'OSi)  =  A£  (13) 

Thus  A£  is  proportional  to  -In  ffip;  the  more  negative 
A£,  the  more  efficient  the  ion  pairing. 

In  principle,  £(M+— OH2)  could  be  obtained  from  the 
mean  value  taking  into  account  the  whole  solvation  shell 
around  the  cation:  -E( =  A H^/n,  where  A 
is  the  enthalpy  of  solvation  and  n  the  coordination  number 
of  the  solvated  cation.  However,  values  of  n  reported  in 
the  literature  depend  strongly  on  the  method  of  mea¬ 
surement.33  Therefore  we  have  estimated  the  energy  of 
the  bond  between  the  cation  and  its  hydrating  waters  from 
partial  charge  calculations  using  a  cation  bound  to  a  single 
water  molecule.  Assuming  an  electron  transfer  involving 
the  cation  M+  and  the  lone  electron  pair  of  the  oxygen 
atom  of  the  water,  the  partial  charge  distribution  and  the 
bond  length  in  the  complex  M+— OH2  were  calculated.  The 
coulombic  term  of  the  bond  energy  was  then  deduced  from 
Ec  *  Q(M)Q(O)/4ir«0r,  where  Q(M)  and  Q(O)  are  the 


(33)  Ami*,  E.  S.;  Hinton,  J.  F.  Solvent  Effect t  on  Chemical  Phenom¬ 
ena;  Academic  Freer  New  York,  1973. 


“All  the  charges  are  given  in  10~3e  (e  *  1.602  x  lCT*  C).  ‘Oxygen 
compiexed  with  the  metal  M.  'Oxygen  in  the  bond  Si-O-Si. 

partial  charges  of  the  cation  and  oxygen,  respectively, 
is  the  vacuum  permittivity,  and  r  is  the  bond  length.  Since 
the  electronegativity  difference  between  the  cation  and  the 
oxygen  atom  is  so  large,  the  Coulombic  contribution  is  the 
only  one  considered. 

The  bond  energy  of  M+-"OSi  was  estimated  by  a  similar 
method.  Using  the  charged  silicate  fragment  shown  in 
Figure  5,  the  cation  is  compiexed  to  the  negative  oxygen 
atoms.  This  structure  occurs  in  basic  aqueous  solutions 
of  silicate  salts  and  consists  entirely  of  Q3  units,  which  are 
probably  similar  to  SiO"  fragments  on  Si02  colloidal  gels.34 
The  hybridization  of  the  oxygen  atoms  has  been  adjusted 
for  a  Si-O-Si  angle  of  about  140°.  The  calculated  partial 
charges  are  shown  in  Table  II. 

The  Coulombic  energies  of  the  bonds  M4"— 'OSi  and 
M+— OH2  are  given  in  Table  III.  The  trend  of  these  en¬ 
ergies  remains  unaffected  either  by  changing  the  number 
of  compiexed  oxygen  atoms  in  the  silicate  fragment  or  by 
using  different  silicate  structures  than  the  anion  consid¬ 
ered.  Since  we  do  not  consider  shielding  and  bonding  by 
the  rest  of  the  solvation  shell,  the  value  of  the  hydration 
energy  is  not  realistic.  However,  the  trend  of  the  hydration 
energy  with  cation  size  remains  meaningful.  The  calcu¬ 
lated  partial  charges  for  the  M+--'OSi  bonds  are  in  rea¬ 
sonable  agreement  with  values  obtained  from  semi- 
empirical  molecular  orbital  methods.36 

The  term  A E  *  -£(M+~'OSi)  +  £(M+~OH2)  of  the 
Gibbs  enthalpy  of  the  ion-pair  reaction  10  was  obtained 
by  using  the  values  for  E(M+— OH2)  and  £(M+— ~OSi)  in 
Table  III  and  is  shown  in  Figure  6.  The  values  of  A £ 
calculated  from  formal  charges  and  from  the  Allred -Ro- 


(34)  Her,  R.  K.  The  Chemistry  of  Silica;  Wiley:  New  York,  1979. 

(35)  KJeinberg,  J.;  Argertinger,  W.  J.;  Griswold.  E.  Inorganic  Chem¬ 
istry;  Heath:  New  York,  1960. 

(36)  De  Jong,  B.  H.  W.  S.;  Brown,  G.  E.  Gtochim.  Coemochim.  Acta 
1*80,  44,  1627. 

(37)  McCormick,  A.  V.;  BeU,  A.  T.;  Radka,  C.  J.  J.  Phys.  Chem.  1M», 
93,  1733. 

(38)  DetelUer,  C.  In  NMR  of  Newly  Accessible  Nuclei;  Unto.  P..  Ed.; 
Academic  Press:  New  York,  1963;  p  105. 
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Table  III.  Estimated  Coulombic  Bond  Energiea  and  Bond  Lesgth*  tor  M*  »  •  •  OHj  and  M*  •  •  •  *OSi 


ionic 


metal 

radius,3 

pm 

£(M*--OH,), 

keal/mol 

r(  M*~OH2), 
pm 

£(M*—'OSi), 

kcal/mol 

pm 

AE. 

kcal/mol 

Li 

60 

53.7 

58.8 

73.7 

184.9 

-20.0 

Na 

95 

27.4 

121.0 

89.5 

204.7 

-42.1 

K 

133 

21.4 

174.4 

66.8 

245.9 

-45.4 

Rb 

148 

20.2 

193.7 

63.7 

265.8 

-43.5 

Cs 

174 

18.3 

222.2 

62.7 

284.6 

-44.4 

*  Fom  Kleinberg  et  al.35 


CMIon  Radius  (pm) 

Figure  $.  Curve  of  A E  vs  ionic  radius  with  A E  calculated  from 
the  bond  orbital  electronegativity  equalization  model. 

chow  model  are  presented  in  Figure  7.  As  can  be  seen  in 
Figure  6,  the  experimental  trend  with  a  maximum  in 
stability  for  potassium  complexation  is  reasonably  well 
predicted  with  the  bond  orbital  electronegativity  equali¬ 
zation  model.  On  the  contary,  if  the  Coulombic  bond 
energies  are  calculated  considering  formal  charges  on  the 
oxygen  atom  and  the  metal,  no  maximum  of  stability  is 
predicted  for  potassium  ion-pair  complexes.  Likewise,  the 
Allred-Rochow  model  does  not  fit  the  experimental  trend 
(Figure  7). 

It  appears  that  the  ion-pair  equilibrium  results  from  the 
balance  between  solvation  of  the  cation  and  strength  of 
the  ion-pair  bond.  For  instance,  though  the  Li+~"OSi  bond 
is  the  strongest,  Li*  is  so  strongly  solvated  that  the  ion 
pairing  is  less  favorable  compared  with  the  other  cations. 
Furthermore,  the  K*— ~OSi  bond  is  more  favored  than 
might  have  been  expected  since  K*  holds  the  hydrating 
water  less  tightly  and  at  the  same  time  the  K*-'OSi  bond 
is  strong  enough.  Potassium  represents  thus  the  best 
compromise  with  its  moderate  size  and  its  low  electro¬ 
negativity.  The  experimental  trend,  and  in  particular  the 
high  stability  of  the  potassium  ion  pair,  is  thus  explained 


Cation  Radius  (pm) 

Figure  7.  Curve  of  AE  vs  ionic  radius  with  AE  calculated  from 
the  Allred-Rochow  model  and  a  formal  charge  model. 

by  our  calculations,  at  least  qualitatively. 

Conclusion 

It  has  been  shown  that  ion-pairing  complexation  be¬ 
tween  cations  and  charged  silicate  groups  occurs  in  basic 
conditions  in  the  condensation  of  TEOS-derived  oligomeric 
structures,  supporting  the  condensation  mechanism  pro¬ 
posed  by  McCormick  et  al.4  NMR  spectroscopy  reveals 
that  the  cation,  initially  solvated  by  water,  can  bind  di¬ 
rectly  to  a  colloidal  silica  particle. 

The  relative  efficiency  of  ion  pairing  and  the  trend  in 
ion-pair  exchange  agree  well  with  a  model  based  on  the 
Bergmann-Hinze  method  of  charge  distribution  calcula¬ 
tions.1  This  reaction  model  assumes  that  charged  silicate 
groups  compete  with  waters  of  hydration  to  coordinate  the 
cation. 
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